HOMOGENIZATION OF A WILSON-COWAN MODEL FOR NEURAL FIELDS 



NILS SVANSTEDT AND JEAN LOUIS WOUKENG 

Abstract. Homogenization of Wilson-Cowan type of nonlocal neural field models is investigated. Mo- 
tivated by the presence of a convolution terms in this type of models, we first prove some general 
convergence results related to convolution sequences. We then apply these results to the homogeniza- 
tion problem of the Wilson-Cowan type model in a general deterministic setting. Key ingredients in this 
study are the notion of algebras with mean value and the related concept of sigma- convergence. 



1. Introduction 

Experiments/observations through different EEG, fMRI, MEG and optical imaging techniques reveal 
electrical activity patterns spanning over several centimeter of brain tissue i.e. of length scales much 
larger than the spatial extent of one single neuron. Moreover, these structures have a lifetime which 
is much larger than the lifetime of typical firing time for a neuron. Now, as the cortex obviously is a 
heterogeneous medium possessing many different spatial and temporal scales it is a need to have rigorous 
ways of determining how the spatio-temporal microstructure is stored in mean field models for the brain 
activity. One way of doing this is by means of homogenization theory based on multi-scale convergence 
techniques. The problem of homogenization or scaling is to determine from data or local characteristics, 
the effective models representing the macroscopic behavior of mesoscopically inhomogeneous media. 

In the quest of studying and understanding the neurodynamics, several heuristic models have been 
developed. The classical leaky integrator unit model [lOl ESI [31] described by ordinary differential equa- 
tions has given rise to the most studied Wilson- Cowan type of models J32] (see also Amari [Tl) which, in 
the one dimensional space, reads as 

d r°° 

—u{x,t) = -u{x,t) + J J{x',x)f{u{x',t))dx'. (1.1) 

Here u{x,t) denotes the neural field which measures the local activity of neurons at position a; G M. The 
integral represents the synaptic input where the function J{x,x') measures the strength of connections 
between neurons at positions x and x' . We refer to this function as the connectivity function. The 
function / is the firing rate function. Equation (|1.1|) models the neural field in a homogeneous medium. 
For more details concerning this model we refer to e.g. Wilson and Cowan [35^ and Amari [I] and the 
more recent work by Coombes [T and the references therein. 

A drawback with the model (|l.ll) above is that does not take into account that the brain is very 
heterogeneous with a structure that exhibits multiple spatial scales ranging from micro- to decimeter. In 
addition the dynamical activities are taking place on multiple time scales. In order to capture several such 
properties we are led to allow the connectivity to depend on both space and time and also to depend on 
multiple spatial and temporal scales. In the present work we will consider time independent connectivity 
kernels. 

An obvious way to impose heterogeneity in the well-accepted Wilson-Cowan or Amari model (jl.ip 
is to replace the connectivity function J by a function J^, where the parameter e > measures the 
heterogeneity of the brain tissue. 
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The homogenization procedure has been used in many apphed science fields to upscale various math- 
ematical models. As far as the neural field models are concerned, there are very few works dealing with 
homogenization techniques. See e.g., [12] or [8]. In [12] a homogenization based approach for study- 
ing nonlocal heterogeneous neural field models of Wilson-Cowan type based on multi-scale expansion 
techniques is initiated. A heterogeneous neural field model is also advocated in [8] where the point of 
departure is the parametrized Wilson-Cowan model 
d f 

— Ue{x,t) = -U,{x,t)+ Je{x - x)f{Ue{x\t))dx\ £ M^, t>Q (1.2) 

Ot J^N 

where Ug denotes the electrical activity level field, / the firing rate function and = Je{x) = J{x,x/e) 
the connectivity kernel which by assumption is periodic in the second argument y — x/ e. In the present 
work, under a general deterministic assumption on the kernel (including the periodicity assumption and 
the almost periodicity assumption) we prove rigorously (see Theorem [5]) that, as e — > 0, the solution 
to the Wilson-Cowan model (|1.2[) converges to the solution ug of a homogenized Wilson-Cowan equation 

d 

g^uo{x, t, y) = -uq{x, t, y) + {J* *f{uoj){x, t, y). (1.3) 
In the special periodic case. Equation (|1.3p reads as 

^Uo{x,t,y) = -UQ{x,t,y)+ [ [ J{x - x' ,y - y')f{uQ{x' ,t,y'))dy'dx' . 

Ot JtsiN Jy 

Due to the nonlinearity in our model equation, we can not use the Laplace transform in the homoge- 
nization process. Also our method works even in the non-Hilbertian framework. Indeed, considering two 
sequences {ue)e and (we)e in i^(M^) and L^iQ) respectively satisfying Ue — >■ uo in L^(K^)-strong S and 
Ve — i> Wo in LP{Q)-weak S as e — J- (where Q is an open subset of M^) we get that Ue * Ve uq * *vq 
in (Q)-weak S as £ — t- 0, where uq * *vo is a double convolution with respect to both macroscopic and 
microscopic variables; see Theorem [5] The above result was first proved by Visintin in the periodic 
setting by using the two-scale transform or unfolding method. Theorem [2] allows us to pass to the limit in 
the convolution terms without using neither the Fourier transform, nor the Laplace transform, and hence 
without restricting ourselves to the Hilbertian setting as it is the case in [3D]. Taking into account the 
fact that the brain is not necessarily a periodic medium (even if it can exhibit some kinds of periodicity) , 
we can therefore emphasize that our work is a true advance in the neural field community. 

The paper is organized as follows. In Section 2 we recall some background material regarding the 
concept of sigma-convergence. We also prove two important results which are of independent interest, a 
general two-scale convergence result for translates (Proposition [3]) and a general two-scale convergence 
result for convolution products (Theorem [2]). The method used in deriving these results is based on 
the notion of algebras with mean value and the concept of sigma-convergence. In Section 3 we prove 
existence of solution to the Wilson-Cowan model and derive the a priori estimate needed for the main 
homogenization result which is stated and proved in Section 4. Finally, Section 5 deals with conclusions 
and outlook. 

2. S-CONVERGENCE AND CONVOLUTION 

2.1. Some properties of algebras with mean value. Let A be an algebra with mean value on 
(see [lTl|34]), that is, A is a closed subalgebra of the C*-algebra of bounded uniformly continuous complex 
functions BUC{R^) which contains the constants, is closed under complex conjugation {u & A whenever 
u € A), is translation invariant (u(- + a) ^ A for any u ^ A and each a £ M^) and such that each element 
possesses a mean value in the following sense: 

(MV) For each u ^ A, the sequence {u^)^yo (where u^{x) = u{x/e), x e R^) weakly *-converges in 
L°°(R^) to some constant function M{u) e C (the complex field). 
It is known that A (endowed with the sup norm topology) is a commutative C*-algebra with identity. 
We denote by A (A) the spectrum of A and by Q the Gelfand transformation on A. We recall that A (A) 
(a subset of the topological dual A' of A) is the set of all nonzero multiplicative linear functionals on 



HOMOGENIZATION OF A NONLOCAL DIFFUSION EQUATION 



3 



A, and Q is the mapping of A into C{A{A)) such that Q{u){s) — {s,u) {s e A{A)), where (,} denotes 
the duahty pairing between A' and A. We endow A{A) with the relative weak* topology on A'. Then 
using the well-known theorem of Stone (see e.g., either [13] or more precisely [H Theorem IV. 6. 18, p. 
274]) one may easily show that the spectrum A{A) is a compact topological space, and the Gelfand 
transformation Q is an isometric isomorphism identifying A with C{A{A)) (the continuous functions on 
A{A)) as C*-algebras. Next, since each element of A possesses a mean value, this induces a mapping 
u I— > M{u) (denoted by M and called the mean value) which is a nonnegative continuous linear functional 
on A with M(l) = 1, and so provides us with a linear nonnegative functional il) i— >■ Mi(?/') — M{Q~^{il})) 
defined on C{A{A)) — Q{A), which is clearly bounded. Therefore, by the Riesz-Markov theorem, Mi(-0) 
is representable by integration with respect to some Radon measure (of total mass 1) in A (A), called 
the M -measure for A 1141. It is evident that we have 



M{u) = / g{u)dp for ue A. (2.1) 
Ja{A) 

The following result is worth recalling. Its proof can be found in [23) , and we recall it here for further 
purposes. 

Theorem 1. Let A be an algebra with mean value on M.^ . The translations T{y) : — > M"'^, T{y)x — 
x + y, extend to a group of homeomorphisms T[y) ; A(A) — > A{A), y G M.^ , which forms a continuous N- 
dimensional dynamical system on A{A) whose invariant probability measure is precisely the M -measure 
13 for A. 

Proof. As A is translation invariant, each translation T{y) induces an isometric isomorphism still denoted 
by T{y), from A onto A, defined by T{y)u = u{- + y) for u e A. Define f{y) : C(A(A)) ^ C(A(A)) by 

T{y)g{u) = g{T{y)u) {u e A) 

where Q denotes the Gelfand transformation on A. Then T{y) is an isometric isomorphism of C{A{A)) 
onto itself; this is easily seen by the fact that Q is an isometric isomorphism of A onto C{A{A)). Therefore, 
by the classical Banach-Stone theorem there exists a unique homeomorphism T{y) of A{A) onto itself. 
The family thus constructed is in fact a continuous A'^-dimensional dynamical system. Indeed the group 
property easily comes from the equality Q{T{y)u){s) = Q{u){T{y)s) {y £ M^, s e A{A), u G A). As far 
as the continuity property is concerned, let {yn)n be a sequence in and {sd)d be a net in A{A) such 
that ?/„ — > y in M.^ and — >■ s in A{A). Then the uniform continuity oi u Q A leads to T{yn)u — T{y)u 
in BUC{R^), and the continuity of Q gives g{T{yn)u) g{T{y)u), the last convergence result being 
uniform in C{A{A)). Hence Q{T{yn)u){sd) Q{T{y)u){s), which is equivalent to G{u){T{yn)sd) -> 
g{u)(T{y)s). As C{A{A)) separates the points of A{A), this yields T{yn)sd T{y)s in A(A), which 
implies that the mapping (y, s) ^ T{y)s, from x A{A) to A{A), is continuous. It remains to check 
that /3 is the invariant measure for T. But this easily comes from the invariance under translations' 
property of the mean value and of the integral representation ()2.1|) . We keep using the notation T{y) for 
T{y), and the proof is complete. □ 

Next, let (1 < P < oo) denote the Besicovitch space associated to A, that is the closure of A with 
respect to the Besicovitch seminorm 

i/p 



limsup-^ / \^iy)fd.y 

r-> + oo |-t>r| J B,. 



where Br is the open ball of M.-^ of radius r. It is known that B^ is a complete seminormed vector space. 
Moreover, we have J5^ C B^ for 1 < p < g < oo. The following properties are worth noticing 124) : 

(1) The Gelfand transformation Q : A C{A{A)) extends by continuity to a unique continuous 
linear mapping, still denoted by Q, of B^ into LP(A(A)), which in turn induces an isometric 
isomorphism Gi, of = B^^/J^ onto LP{A{A)) (where JV ^ {u e B^ : g{u) = 0}). Furthermore 
if u e n L°°(R^) then g{u) e L°°(A(A)) and ||e(u)||ioc(A(A)) < \Ml^(r«)- 
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(2) The mean value M viewed as defined on A, extends by continuity to a positive continuous linear 
form (still denoted by M) on satisfying M{u) ~ J^^^.^Q{u)dl3 {u e B^). Furthermore, 

M{u{- + a)) = M{u) for each u S and aU a e M^. 

(3) The dynamical system T(y) is ergodic if and only if for every u G B\ such that \\u — u{- + y)\\p — 
for every y e we have ||m — M{u)\\p — 0. 

In order to simplify the text, we will henceforth use the same letter u (if there is no danger of confusion) 
to denote the equivalence class of an element u g -B^. The symbol g will denote the canonical mapping 
of onto = B^/Af. For u G (resp. u e B^) we shall set u = Gi{u) (resp. u = G{u))- 

In our study we will deal with a special class A of algebras with mean value. An algebra with mean 
value A is in A if it is separable (hence its spectrum A{A) is a compact metric space) and further A{A) has 
group property. Here below are some examples of algebras with mean value verifying the above property. 
1) Let A — Cpci{Y) {Y — (0, 1)^) denote the algebra of F-periodic complex continuous functions on 
M.^ . It is well-known that A{A) = R^/Z^ (the iV-torus) which is a separable compact metrizable 
topological group. 2) If we denote by AP(R^) the space of complex continuous almost periodic functions 
on [1[S], then for any countable subgroup U of R^, denoting by APk(R^) the subspace of AP(R^) 
consisting of functions that are uniformly approximated by finite linear combinations of functions in the 
set {7^ : k e TZ} (where 7^.(2/) = exp(2i7rfc ■ y), y e R^), we get that APtz{R^) e A with the further 
property that its spectrum is a compact topological group homeomorphic to the dual group TZ oi TZ 
consisting of characters 7^. (fc 6 TV) of R^; see [T5j for details. 3) Finally, let i3oo(R^) denote the space 
of continuous complex functions on R^ that have finite limit at infinity. It can be easily shown that its 
spectrum consists of only one point, the Dirac mass at infinity (5 00, hence is a compact topological group. 
Also Soo(R^) is separable (see [2]), so that Soo(R^) G A. Some other examples can be considered by 
taking a combination of the previous ones. 

Some notations are in order. Let A G A. Since A(A) is a topological group (which we henceforth 
denote additively), the mapping (s,r) s + r is continuous from A(A) x A(A) into A(j4). ~s shall 
stand for the symmetrization of s £ A(A). Now, in the above notation, if we take a look at the proof of 
Theorem [T] we observe that the dynamical system associated to the translations is defined by 

T{y)s ^5y + s for (y, s) G x A(A) 

where 5y is the Dirac mass at y. With this in mind, for s — 5x and r — Sy we may see that we have in 
the same notations, s + r = dx+y and s — r ~ Sx-y 
We get the following result. 

Proposition 1. Let A e A. Let u e A and a G R^. Then 

giu{-+a)) ^g{u){- + 5a). (2.2) 

Proof. Let us recall that (|2.2p is equivalent to G{u{- + a)){s) — g{u){s + 6a) for any s G A{A). So, let 
y eR^, then 

g{u{- + a)){6y) = {Sy,u{- + a)) = u{y + a) 
= {Sy+a,u) = Q{u){dy+a) 

= g{u){6y+6a). 

Now, let s G A(A); there exists a sequence (j/„)„ C R^ such that Sy^ ^> s in A{A) as n ^ 00 (indeed 
the set {6y : y G M^} is dense in A{A) which is a metric space). Since A(A) is a topological group, 
the mapping (s,r) ^ s + r is continuous from A{A) x A{A) into A(A). So for fixed r — 6a we have 
that Sy^ + 5a — > s + 5a, and by the continuity of G{u), it comes that Q{u){5y^ + 5a) — > G{u){s + 5a) 
and g{u{- + a)){5y^) — ;> g{u{- + a)){s) as n — >■ 00. The uniqueness of the limit gives by the fact that 
Giui- + a))i5yj =g{u){5y,^ + 5a), GH- + a))(s) = e(u)(s + 5a). □ 
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2.2. E- convergence and convolution results. We begin with the definition of the concept of S- 
convergence. Let A € A, and let Q be an open subset of R^. Throughout the paper the letter E will 
denote any ordinary sequence E = (£„) (integers n > 0) with < e„ < 1 and £„ — as n — ^ oo. 

Definition 1. (1) A sequence (we)^>Q C (Q) {1 < p < oo) is said to weakly 'S-converge in (Q) to 
some uq G LP{Q; B^) if as e — 0, 

/ Me {x) tp^ {x) — ?• / / Uq {x, s) ip (x, s) dxdfH (s) (2-3) 

Jq JJqxA{A) 

for all ■0 G i'' (Q; A) {1/p' = 1 — 1/p) where tp^ (x) = tp (x, x/e) and ^p (x, •) — G{ip (x, •)) a.e. in x £ Q. 
We denote this by — J> uq in L'°(Q)-weak S. 

(2) A sequence (ue)e>o C LP{Q) (1 < p < oo) is said to strongly Yi-converge in LP{Q) to some uq G 
L^{Q] /B^) if it is weakly E-convergent and further satisfies the following condition: 

\We\\LP{Q] II"oIIlp(QxA(A)) • (2-4) 

We denote this by — J> in LP((5)-strong S. 

We have the following result whose proof can be found in [IB] (see also [5] ) . 

Proposition 2. (i) Any hounded sequence {u^)^^e in LP{Q) (where 1 < p < oo) admits a subsequence 
which is weakly Yi-convergent in LP{Q). 

(ii) Any uniformly integrable sequence {u^)^(ze in L^{Q) admits a subsequence which is weakly Yi-convergent 
in L\Q). 

We recall that A sequence (w£)e>o in L^iQ) is said to be uniformly integrable if it is bounded in L^{Q) 
and further satisfies the property that supgj,g \ue\dx — > for any integrable set X C Q for which 
\X\ — > 0, where \X\ denotes the Lebesgue measure of X. 

Now, fix t € M^. Then {St/e)e>o is a sequence in the compact metric space A(yl), hence it possesses 
a convergent subsequence still denoted by {St/e)s>o- In the sequel we shall consider such a subsequence. 
Let r G A (A) be such that 

(5i -> r in A{A) as e ^ 0. (2.5) 

Finally let Q be an open subset in M^, and let {u^)^^o be a sequence in LP{Q) {1 < p < oo) which is 
weakly S-convergent to uq G LP{Q; B^). Define the sequence {ve)e>o as follows: 

Vs{x) = u^{x + t), x £ Q — t. 

We then get the following result. 

Proposition 3. Assume (|2.5p holds. Then as e — > 

Ve vq in U\Q — t)-weak E 
where vq G U'{Q ~ t, S^) is defined by vq{x, s) = uq{x + t,s + r) for (x, s) G (Q — i) x A(yl). 
Proof. Let if G C^{Q — t) and ip £ A. Let {yn)n be an ordinary sequence (independent of e) such that 

5y^ — r in A(^) as n ^ oo. 
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We have 



Q-t 

Ue{x)(p{x - t)lp 
t) 



Q 



= / Ue{x)ip{x 
JQ 



£ 

X t 

e e 
X t 
e £ 



dx 



>(f-yn) 



dx 



+ j Ue{x)ip{x - t)tp (^^ - dx 



= {i) + {n). 



On one hand, as s — >^ 0, 



(//) -> // uq{x, s)ip{x - t)g{iP{- - yn)){s)dxdl3{s). 

J JQxA(A) 



/QxA(A) 

But we have g{tp{- — Vn)) = ~ ^yn) V'(' ~ ^) uniformly in C(A(^)) as n — >■ oo. Hence 



QxA(A) 



Uo{x, s)if{x - t)g{ip{- - yn)){s)dxdP{s) 



QxA(A) 



uo{x, s)(fi{x — t)ip{s — r)dxdj3{s) as n — >■ cx), 



and 



uq{x, s)(p{x — t)ip{s — r)dxdj3 



'QxA(A) 

On the other hand, 



uo{x + t,s + r)(p{x)ip{s)dxdj3 . 



(Q-t)xA(A) 



< c sup 



\Ue{x)\ lifiix-t)] 
t 



X t 

e e 



i) z 



tp{z- Vn) 



dx 



Vn) 



Now, using the uniform continuity of ip, we obtain 
It therefore follows that, as £ ^ 0, 



{Q is an isometry) 



as £ — >■ and next n — > oo. 



IQ-t 

This concludes the proof. 



/ Ue{x + t)ip{x)ip (—j dx ^ // uo{x + t,s + r)(p{x)il){s)dxd^. 

JQ-t J J(Q-i)xA{A) 



□ 



The next important result deals with the convergence of convolution sequences. Let p > 1 be a real 
number, and let {ue)e>Q C LP{Q) (where we assume here Q to be bounded) and (i'e)e>o C L-^(M^) be 
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two sequences. One may view as defined in the whole by taking its zero extension off Q. Define 

{u,*v,){x)^ / u,{t)v,{x -t)dt (a;6R^). 
For u G LP{M.^; B^) and v G L^{M.^ ; B\) we define the convolution product u* *v as follows 
Qi{u * *v){x, s) := // u{t,r)v{x — t, s — r)dtdl3{r) 

J JtS."xA(A) 

for {x, s) G X A{A). We denote Qi {u * *v) by {u * *v). This gives a function in the space LP(M.^ ; B^) 
with the property 

\\U * *f llLP(RiV.gP) < l|M||ip(aN.gP ) ||w||2,l(RiV.gl^) • 

The above inequality can be checked exactly as the Young inequality for convolution. We have the 
following result. 

Theorem 2. Let (u^)e and (we)^ be as above. Assume that, as e — s- 0, — J- uq in LP{Q)-weak E and 
We Vq in L^{M.^) -strong S, where uq and vq are in LP{Q;B^) and L^{M.^;B]^) respectively. Then, as 

Me * We — > Uq * *va in L^{Q)-weak S. 
Proof. Let 77 > and let ipQ G IC{M.^;A) (the space of continuous functions from M.^ into A with 



compact support) be such that 



Li(R"xA(A)) 



< ^. Since vq in L-'^(R )-strong S we have 



Wo - -00 

that ~ i/jq Vq — i/jq in i-'-(K^)-strong S, hence ||we — ■'/'oIIli(r") ~^ ^0 — V'o 
So, there is a > such that 

ll^'e - 0oIIli(b«) < ?7 for < £ < a. 
For / G IC{Q; A), we have (by still denoting by the zero extension of off Q) 



Li(R"xA(A)) 



as e — > 0. 



(2.6) 



(ue * v^){x)f {x, da; 



Ue{t) 
Ue{t) 
Ue{t) 
Ue{t) 

/r« 

{I) + ill)- 



:{t)v,{x~t)dtj f (^X,^'^ 
Veix - t)f (^X, ^) 



£ 

dx 



dx 
dt 



Ve{x)f [ X + t, — ^ — ] dx 



(veix) ~ ii%{x))f{x + t)dx 
il)%{x)f{x^t)dx \ dt 

R" / 



dt 
dt 



On one hand one has (/) = JqIu^ * (we — 4'o)]{x)f^{x)dx and 

l(/)| < hellLP(Q) - V-oIIli 
< C ||Ue — V'oIIl1(R«) 

where c is a positive constant independent of e. It follows that 

(/)| < C77 for < £ < q; 
On the other hand, in view of Proposition [3l we have, as £ ^> 0, 



\LP'(Q) 



(2.7) 



^Pl{x)f{x + t)dx 



i^l{x~t)f{x)dx 

iPq{x — t,s — r)f{x, s)dxd/3{s), 

R"xA(A) 
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where r — \m\5t/e (fo^' '"^ suitable subsequence of £ — 0) in A(yl). So let $ 

by 

r) = 1 1 ijo{x-t,s- r)f{x, s)dxdP{s), {t, r) G x A(A). 

J JR"xA(yl) 

Then it can be easily checked that $ G /C(R^; C(A(A))), so that there is a function ^E* G 
$ = C/ o We can therefore define the trace "^^{t) ^(t^t/e) [t G M^) and get 

^%t) = (5,,e(*(t,-))) 



A(A) ^ M be defined 



; A) with 



V'ola^ — s — S±) f{x, s)dxdl3{s). 



Next, we have 
(.II) 

As for {III), set 



"xA(yl) 



Ue{t){ i^l{x)f{x + t)dx^'$''{t)]dt+ Ue{t)-^^{t)dt 



(//l) + (7/2). 



K W = / ^l{x)r-{x + - for a.e. t G 



We claim that, for a.e. t, Vi(i) — J- as e — J- (possibly up to a subsequence). Indeed, due to Proposition 

4'q{x) {x + t)dx jj 'iIjq{x — t,s — r)f{x,s)dxdp{s) as e ^ 

J JK"xA(A) 

where r is such that 6t/^ r in A{A) for some subsequence of e. Moreover, since ^''^(i) = ^'^{t, 5t/e), by 
the continuity of •), we have for the same subsequence as e ^> that 



■0o(a^ — t,s ~ r)f{x, s)dxdp{s). 



'xA{A) 

The above claim is justified. Furthermore 

\Ve{t)\ < c for a.e. t G Q 

where c is a positive constant independent of t and e. On the other hand, since Q is bounded, we infer 
from the weak S-convergence of ('Ue)^ that ~^ Ja{A) ^o('' T)dl3(r) in L^(Q)-weak as e — ?> 0. Therefore, 
by fS", Lemma 3.4] it follows that (J/i) as e 0. 

Regarding {II2), using once again the weak S-convergence of (ug)^, we get 



uo{t,r)'^{t,r)dtdl3, 



QxA{A) 



and 



uo{t,r)'i!{t,r)dtdl3 



QxA{A) 



QxA{A) 



uo{t,r)<^>{t,r)dtdl3 



ua{t, r)ipQ{x — t,s — r)dtdl3(r) 



QxA(A) \ J JR" xA(yl) 

(wo * *V'o)(2;, s)f{x, s)dxd(i{s) 

QxA(A) 



/(x, s)dxdl3{s) 
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Thus, there is < ai < a such that 

Q JJqxA{A) 
Now, let < e < ai be fixed. From the following decomposition 



< ^ for < £ < ai. 



(2.8) 



(Ue * Ve)f^dx - 

)xA(A) 
[u^ * (We - V'o)] fdx + 



+ / (We * i^Dfdx 

JQ 

we infer from ([^ -(1^ that 

(lie * Ve)rdx - jj 



m * ^itpQ - vo) 

QxA{A) L 

{uq * *ipQ)fdxdl3, 



fdxd/3 



QxA{A) 



{uq * *vo)fdxd/3 



Q 



xA{A) 



< Crj for < £ < ai 



Here C is a positive constant independent of e. This concludes the proof. 



□ 



Remark 1. In this work, we will deal with sequences of functions (ite)e>o in the space LP{Qt) (where 
Qt = Q X (0,T)), and we will say that such a sequence is weakly S-convergent in LP{Qx) if 



{x,t) f (x,t,^'j 



dxdt 



QtxA(A) 



Uq (x, t, s) / (x, t, s) dxdtd/S (s) 



for all f e LP (Qt;^) (l/p' = 1 - where fix,t,-) = g{f{x,t,-)) a.e. in {x,t) e Qt- With the 
above definition, it is a very easy exercise to see that all the previous results of the current subsection 
are carried over mutatis mutandis to the present setting. 

3. Statement of the problem: existence result and a priori estimate 
We consider the parametrized Wilson-Cowan model [35] 



ite(a:, 0) — M°(x), X € . 



,ix, t) + 4„ J%x - 0/ [l,u,i^, i)) d^, x e M^, t > 

^ ^ r- TON 



(3.1) 



where ite denotes the electrical activity level field, / the firing rate function and J'^ = J'^{x) = J{x,x/e) 
the connectivity kernel. We assume that J S /C(M^;^) (where ^ G A) is nonnegative and is such that 



/rn J^{x)dx < 1, / : X M is a nonnegative Caratheodory function constrained as follows: 

(HI) For almost all y G M^, the function /(i/, ■) : A n- /(y. A) is continuous; for all A € R, the function 
/(•, A) : y f{y,X) is measurable and f{-,0) lies in L^(R^) n L^(R^); there exists a positive 
constant ki such that 

\f{y,fJ-i) - fiy,f^2)\ < fci IMi for all y e and aU jii,ji2^ ^■ 

An example of a function / that satisfies hypothesis (HI) is f{y,X) — g{y)h{\) where g £ /C(R^), 
g > and 



h{X) 



(A e R) where /3 > and 6* are given. 



l + exp(-/3(A-6i)) 

It follows from (HI) that, for any function u £ L^(R^), the function cc 1— t- f{x/e,u{x)) denoted below by 
/^(•, It), is well defined from R^ to L-^(R^). Moreover it is an easy exercise to see that 

||r(-,M)|li2(R«) < fci ||w||^2(R«) +ci forallO<£<l (3.2) 

where ci = ll/(': 0)IIl2(rjv). The following existence result holds true. 
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Theorem 3. Let < T < oo. Assume that e L^{M.^) n L'^{K'^). Then for each fixed e > 0, there 
exists a unique solution G C([0, oo); i"'^(R^) O.L'^iM.^)) to (I3.ip which satisfies the estimate 

We(-.OIlLi(R«) + ||we(-,0llL2(R") <Cforalle>0 (3.3) 



sup 

0<t<T 

where C is a positive constant depending only on u'^ and T . 

Proof. Inspired by the proof of Theorem 2.1 in [17] (see also [25]) we define the space 

X = C([0,p];Li(R^)nL2(MJV)) 

with norm 

\W\\x = max (||w(-,t)||^i(]g«) + ||u(-,i)llL2(E")) : 

where p > is an arbitrary real number to be suitably chosen in the sequel. Next, by the assumptions 
on J we have that the trace is well defined and is an element of /C(K^) for every e > 0. Let us now 
consider the operator : X ^ X defined by 

KM{x,t)^4>{x,Q)+ I {J' *f{;cj,)-cj)){x,T)dT. 



We observe that a fixed point to i^e is a local solution £ X to p.ip by the Banach fixed point theorem. 
We proceed by showing that is a strict contraction on X. By the properties of J^, / and by the use 
of Young's inequality, we obtain 

||A',(u)-i^,(«)||^.(R„)(t) 



< 



* (/"(•:") -/^(•,«)IIl2(js.«) (t) + w||i2(R«) {T)]dT 



P 



We also have 



< / (fcl + l)||lt-v||^2(RiV)(T)dT 
•J 

< (fci + l)p \\u — v\\^ for every < t < p. 

\Ke{u) - Ke{v)\\^i,jgN) (t) < (fci + l)p\\u - v\\x for every < t < p. 



Choosing p small enough such that 2(fci + l)p < 1, we obtain that is a strict contraction on X. Thus 
by the Banach fixed point theorem, there exists a unique local solution £ X to (jSTJ). Next, arguing 
exactly as in the proof of Theorem 2.7 in [19], we get the global existence of the solution to (|3.ip . Now, 
we need to check that estimate (|3.3p holds true. In order to do that, by multiplying Eq. (jSHJ) by Ue{x, t) 
and integrate the resulting equality over R^, we obtain 

^/ \u^{x,t)\^ dx ^ -2 [ \u,{x,t)fdx + 2 [ Ueix,t)[r * f%-,Ue)]ix,t)dx. 
Ot Jrn Jrn Jrn 

By Holder and Young inequalities, we get 

/ Ueix,t)[r *f{;Ue)]ix,t)dx < (', || i2(R«) ( / * f {■ , U,)\'' dx] 

^ 0IIl2(jjJV') ||>^^|lil(RiV') ||/^('lUe)lli2(RjV') . 

We infer from ((3?2]) that 

lir(-:We)||i2(RiV) < fci ||Me(-,i)llL2(R«) +Ci. 
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Thus 



^ II^e('; Oil L2(Riv-) < —2\\u^{-,t)\\ 2;,2(KN) + 2fci i)|| £,2(RJV) + 2ci i)|| i2(RjV) 



2(R") + ||w£(-,0IIl2(R«) 

Either we have 



II"£(-'0IIl2(rjv) < — for all e > 



2ci 

or there are some e > for which 

II / Ml ^'^1 
lFe(')*)llL2(RN-) > — . 



\\Ue{-,t)\\l2(jiN) < 2 ||We(-,i)lli2(RiV) (^fcl + = 3fcl ||m£(- 



|2 



Then, for such e we have 
hence 

I|w£(-:0IIl2(rjv) < exp(3fcii) ||u°||^2(rjv) < exp(3fcir||u°||^2(Riv) • 
In both cases we see that there exists a positive constant C depending only on both u° and T such that 

||ue(-,i)|lL2(RiV) < C for all < e < 1 and all < t < T. 
It is also an easy task to see that 

\We{-,t)\\^i^^N^ < C for all e > and all < t < T. 
This completes the proof. □ 

Remark 2. From the uniform boundedness of (ue)o<e<i in C([0, T];L^{R^) D L'^{R^)), we deduce that 
(ue)o<e<i is uniformly integrable in L^{R^ x (0, T)). Indeed, let B C x (0, T) be an integrable subset. 
Denoting by \B\ its Lebesgue measure, we have by Holder's inequality that 



/ \u,\dxdt < \B\^\\u,\\^, 

JB 



2(R«x(0,T)) 

< C\B\^, 

hence supQ<g<]^/g \ue\dxdt — > when \B\ — > 0. Wc may therefore use [part (ii) of] Proposition [2] to 
deduce the existence of a subsequence of (ue)o<£<i that weakly E-converges in L^(M^ x (0,r)). 

4. HOMOGENIZATION RESULT 

Let A be an algebra with mean value taken in the class A. In order to perform the homogcnization 
process, we assume that the function / satisfies the following hypotheses (in which we set = x 
(0,T)): 

(H2) f{-,n) g A for all e K 

(H3) For any sequence {ve)e>o C ^-'^(IR^) such that — > vq in L-'^(M^)-weak S, wc have f^{-,Ve) ~^ 
.f{-,vo) in Li(R^)-weak E. 

The hypothesis (H3) is meaningful. Indeed the convergence result vq in L^(M^)-weak S does 

not entail the convergence result f{-,vo) in L^(M^)-weak S in general. However there are 

many situations in which assumption (H3) is satisfied. Here below are some examples. 

1) Assume that A = Cpci{Y) and that f{-,n) is F-periodic, f{y, •) is convex. Assume further that 
/ satisfies: 
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(H4) f {j,Ve{x,t)) dxdt JJ^N fiy^'"aix,t,y))dxdtdy as e 0, whenever Vq in 

LH]R^)-weak S. 

Then /''(•, We) fi-,vo) in L^{R^)-weak E as £ -> 0; see [H Theorem 4.3]. 
2) Assume hypotheses (H1)-(H2) hold true. Let (we)e>o be a bounded sequence in L^{M.^)nL'^{M.^) 
such that Ve — wo in L'^{M.^)-weak S. Then (possibly up to a subsequence) we have —^■ 
vq in L^(R.^)-weak E. In view of assumption (HI), the sequence {f^ {■,v^))^^q is bounded in 
L-'^(R^) n L^(R^) so that, up to a subsequence we have /^(•, Wg) zq in L^(M^)-weak E (where 
zq G i^(M:^; ;6^)). Now, if we further assume that 



liminf / {■,Ve)vedxdt < If ZQVQdxdtdfi, 



;xA(A) 

then f''{-,Ve) -)> /(-.i^o) in Li(M^)-weak E; see [33l Theorem 8]. 
3) In the special case when /(y, A) = g{y)h{\) with g G An/C(M^), 5 > and h{\) = ij^cyi-p(~p(\~e)) ' 
(H3) is still true. Indeed, as stressed in Appendix A of j5], we have h{v^) -4 ft.(wo) in L^(M^)-weak 
E, so that, since g £ A, /^(•,Wg) — g^h{v^) — > gh{vQ) = /(•,wo) in i"'^(MiJ^)-weak E. 
We can now state and prove the homogenization result. 

Theorem 4. For any fixed e > 0, let he the unique solution of p.ip . Then as e — > 0, we have 

Ue uo in L^{R^)-weak E (4.1) 

where uq G C{[0,T]; L^(R^ ; B\)) is the unique solution to the following equation 

^{x, t, y) = -Mo(x, t, y) + {J* */(•, uo)){x, t, y), (x, t)&R^, yeR^ 
uo{x, 0, y) = u°{x), x&R^, y G E^. ^ ' 

Proof. We infer from Remark [5] that the sequence (ue)e>o is uniformly integrable in L^{R^). So, given 
an ordinary sequence E it follows from [part (ii) of] Proposition [5] that there exist a subsequence E' of 
E and a function uq G L^(R|(;i3]4) such that, as i?' 3 e — >■ 0, we have (|4.ip . It now remains to check 
that itQ solves Eq. (|4.2p . Indeed it can be easily shown that in view of the properties of / and J, the 
solution to (|4.2p is unique, so that, by the uniqueness property we have the convergence result ()4.ip for 
any ordinary sequence E (not only up to a subsequence E'), and hence for the whole sequence e ^ 0. 
Now, since J G IC{R^ ; A) C /C(M^ x (0, T); ^), we have that 

^ J in L^{R^ X (O,r))-strong E as -B' 9 £ 0; see e.g. [H]. 

Hypothesis (H3) together with the convergence result (14. ip lead us to 

/^(•,We) ^ /(^wo) in L\R^)-weak T. as E' B e ^ 0. 

It therefore follows from Theorem [5] and Remark [T] that 

J" * f'{-,ue) -> J* *f{-,uo) in Li(R^)-weak E as £" 3 £ ^ 0. 

Next, Eq. p.ip is equivalent to the following integral equation 

u,{x, t) = «°(a;) + / [(r- * r{; uM^, t) - u,(x, r)] dr. 
Jo 

Hence, letting E' 3 e ^ Q and using Fubini and Lebesgue dominated convergence results in the integral 
term, we end up with 

uo(a;,i,2/) = M°(a;) + / [(J **/(•, Mo))(a;, r, y) - uo(a;, r, y)] dr 
Jo 

which is equivalent to (14. 2p . Moreover this shows that mq lies in C([0, T]; i^(R^; B]^)) as expected. This 
concludes the proof. □ 
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5. Conclusions and outlook 

In this paper we have proved some important results which are relevant to the theory of homogenization 
in connection with convolution sequences (see e.g. Theorem[2]). It is to be noted that Theorem [2] general- 
izes to the case of algebras with mean value, its counterpart proved by Visintin [29] in the special context 
of the algebras of continuous periodic functions. This result, based on the so-called sigma convergence 
concept, has allowed us to efficiently upscale a heterogeneous Wilson-Cowan type of models for neural 
fields. The homogenization result derived is accurate and more likely can not be achieved through other 
classical and conventional methods such as the asymptotic expansions or the time averaging. Another 
aspect that should be emphasized is that Theorem [2] widely opens the door to many other applications 
in applied science. In this regard, it may allow to study homogenization problems in connection with 
partial differential equations with fractional order derivatives. Such kinds of homogenization problems 
have not yet been solved so far. 

Acknowledgement 1. The authors would like to thank the anonymous referee for valuable remarks 
and suggestions that helped them to significantly improve they work. 
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